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We consider phase turbulent regimes with nonzero winding number in the one-dimensional Com¬ 
plex Ginzburg-Landau equation. We find that phase turbulent states with winding number larger 
than a critical one are only transients and decay to states within a range of allowed winding numbers. 
The analogy with the Eckhaus instability for non-turbulent waves is stressed. The transition from 
phase to defect turbulence is interpreted as an ergodicity breaking transition which occurs when the 
range of allowed winding numbers vanishes. We explain the states reached at long times in terms 
of three basic states, namely quasiperiodic states, frozen turbulence states, and riding turbulence 
states. Justification and some insight into them is obtained from an analysis of a phase equation for 
nonzero winding number: rigidly moving solutions of this equation, which correspond to quasiperi¬ 
odic and frozen turbulence states, are understood in terms of periodic and chaotic solutions of an 
associated system of ordinary differential equations. A short report of some of our results has been 
published in [Montagne et ah, Phys. Rev. Lett. 77, 267 (1996)]. 

PACS: 05.45.-fb,82.40.Bj,05.70.Ln 


I. INTRODUCTION 


A. The complex Ginzburg-Landau equation and its phase diagram 


Spatio-temporal complex dynamics is one of the present focus of research in nonlinear phenomena. This 

subject lies at the intersection of two important lines of thought: on the one hand the generalization of the ideas of 
dynamical systems theory to high dimensional situations and on the other the application of some concepts 

and tools developed in the field of statistical mechanics, specially in the study of phase transitions, to the analysis of 
complex nonequilibrium systems 0-§- 

An important effort has been devoted to the characterization of different dynamical states and transitions among 
them for model equations such as the Complex Ginzburg-Landau Equation (CGLE) 
equation for a complex field A(x, t): 


10- im. The CGLE is an 


dtA = A -I- (1 -I- ici)V^A — (1 -b ic2) | A p A . 


( 1 . 1 ) 


A(x, t) represents the slowly varying, in space and time, complex amplitude of the Fourier mode of zero wavenumber 
when it has become unstable through a Hopf bifurcation (the signs used in dO) assume it to be supercritical). 
The CGLE is obtained univ ersally when analyzing the dynamics sufficiently close to the bifurcation point. In one 
dimensional geometries, ( |l.l| ) or a coupled set of similar equations with additional group velocity terms describe also 
the evolution of the amplitudes of Hopf-bifurcated traveling waves p|, p^ , ^ . Binary fluid convection , transversally 
extended lasers po| , ^ , chemical turbulence bluff body wakes the motion of bars in the bed of rivers 

[p^ , and many other systems have been described by the CGLE in the appropriate parameter range. We will restrict 
ourselves in this paper to the one-dimensional case, that is A = A(x, t), with x G [0, L]. As usual, we will use periodic 
boundary conditions in x. 

The one-dimensional Eq. (0) has traveling wave (TW) solutions 

Afe = v'l - Uk = C 2 A- (ci - C 2 )fc" (1.2) 

with k G [—1,1]- When 1 -b ciC 2 > 0 there is a range of wavenumbers [—fc_E,A:_E] such that TW solutions with 
wavenumber in this range are linearly stable. Waves with k outside this range display a sideband instability (the 
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Eckhaus instability |l| Jl3| , |26| ). The limit of this range, vanishes as the quantity 1 + C 1 C 2 approaches zero, so that 
the r ang e of stable traveling waves vanishes by decreasing 1 + C 1 C 2 . The line 1 + C 1 C 2 = 0, is the Benjamin-Feir-Newell 
line |27y2^ , labeled BEN in Fig. 0 Above that line, where 1 + C 1 C 2 < 0, no traveling wave is stable and different 
turbulent states exist. A major step towards the analysis of phases and phase transitions in (1.1) was the numerical 
construction in [fi], pd]jl^ ] of a phase diagram that shows which type of regular or chaotic behavior occurs in different 
regions of the parameter space [ci,C 2 ]. Fig. has been constructed from the data in Above the BEN line, 

three types of turbulent behavior are found, namely phase turbulence (PT), defect or amplitude turbulence (DT), 
and bichaos (BC). 



Cg 


FIG. 1. Regions of the parameter space [ci,C2] for the CGLE displaying different kinds of regular and chaotic behavior. 
Lines Li, L3 were determined in ||ij, |ll| , [l2| . See the text in Section II for the explanation of the different symbols. 

Phase turbulence is a state in which A{x,t) = |A|e*'^ evolves irregularly but with its modulus always far from 
|A| = 0. Since the modulus never vanishes, periodic boundary conditions enforce the winding number defined as 

iy = ^ j da:(pdx (1.3) 

to be a constant of motion, fixed by the initial condition, n is always an integer because of periodic boundary 
conditions . The quantity k = 2 'ki' jL can be thought of as an average or global wavenumber. To the left of line Li 
(region DT), in contrast, the modulus of A becomes zero at some instants and places (called defects or phase slips). 
In such places the phase ip becomes undefined, so allowing 1/ to change its value during evolution. BC is a region in 
which either PT, DT, or spatial coexistence of both can be observed depending on initial conditions. It should be 
noted that chaotic states exist also below the BEN line: To the left of the line L2, a chaotic attractor called Spatio 
Temporal Intermittency (STI) coexists with the stable traveling waves |j^. A diagram qualitatively similar to Fig. 
has also been found for the two-dimensional CGLE |^J3^ . Despite the relevance of v in the dynamics of the CGLE, 
most studies of the PT regime have only considered in detail the case of n = 0. In fact the phase diagram in Fi g. ^ 
was constructed 000 using initial conditions that enforce i/ = 0. Apart from some limited observations |0 , 0 , pO| ( 
systematic consideration of the yf 0 {wound) disordered phases has started only recently . States with 

n ^ 0 are precisely the subject of the present paper. 


B. The PT-DT transition 

Among the regimes described above, the transition between PT and DT has received special attention 
[f7|, pT|Jl^JM| -^ . The PT regime is robustly observed for the large but finite sizes and for the long but finite ob¬ 
servation times allowed by computer simulation, with the transition to DT appearing at a quite well defined line {Li 
in Fig. 0 ) 00, but it is unknown if the PT state would persist in the thermodynamic limit L —> 00 . One possible 
scenario is that in a system large enough, and after waiting enough time, a defect would appear somewhere, making 
thus the conservation of v only an approximate rule. In this scenario, a PT state is a long lived metastable state. In 
the alternative scenario, the one in which PT and the transition to DT persist even in the thermodynamic limit, this 
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transition would be a kind of ergodicity breaking transition in which the system restricts its dynamics to the 

small portion of configuration space characterized by a particular v. DT would correspond to a “disordered” phase 
and different “ordered” phases in the PT region would be classihed by its value of v. The idea of using a quantity 
related to v as an order parameter @ has also been independently proposed in |3lj. 

The question of which of the scenarios above is the appropriate one is not yet settled. Recent investigations seem 
to slightly favor the first possibility [^p'slJl^JdOtl . The most powerful method in equilibrium statistical mechanics to 
distinguish true phase transitions from sharp crossovers is the careful analysis of finite-size effects Such type of 
analysis has been carried out in |l5y3C| ] , giving some evidence (although not definitive) that the PT state will not 
properly exist in an infinite system or, equivalently, that the Li line in Fig. approaches the BFN line as L —> oo. 
Here we present another finite-size scaling analysis, preliminarily commented in [|l0| , based on the quantity v as an 
order parameter. Our result is inconclusive, perhaps slightly favoring the vanishing of PT at large system sizes. In 
any case, the PT regime is clearly observed in the largest systems considered and its characterization is of relevance 
for experimental systems, that are always hnite. In this paper we characterize this PT regime in a finite system as 
we now outline. 


C. Outline of the paper 


We show that in the PT regime there is an instability such that a conservation law for the winding number occurs 
only for v within a hnite range that depends on the point in parameter space. PT states with too large \v\ are only 
transients and decay to states within a band of allowed winding numbers. Our results, presented in Section allow 
a characterization of the transition from PT to DT in terms of the range of conserved v. as one moves in parameter 
space, within the PT regime and towards the DT regime, this range becomes smaller. The transition is identihed 
with the line in parameter space at which such stable range vanishes. Analogies with known aspects of the Eckhaus 
and the Benjamin-Feir instabiliti es a re stressed. States with v ^ Q found in the PT region of parameters at late times 
are of several types, and Section HI describes them in terms of three elementary wound states. Section fV gives 
some insight into the states numerically obtained by explaining them in terms of solutions of a phase equation. In 
addition, theoretical predictions are made for such states. The paper is closed with some final remarks. An Appendix 
explains our numerical method. 


II. THE WINDING NUMBER INSTABILITY 


The dynamics of states with non-zero winding number and periodic boundary conditions has been studied numer¬ 
ically in the PT region of parameters. In order to do so we have performed numerical integrations of Eq. (1.1) in 
a number of points, shown in Fig. |^. Points marked as O correspond to parameter values where intensive statistics 
has been performed. The points overmarked with x correspond to places where finite-size scaling was analyzed. 
Finally the symbol -I- correspond to runs made in order to determine accurately the PT-DT transition line (Li). Our 
pseudospectral integration method is described in the Appendix. Unless otherwise stated, system size is L = 512 and 
the spatial resolution is typically 512 modes, with some runs performed with up to 4096 modes to conhrm the results. 
The initial condition is a traveling wave, with a desired initial winding number i>i, slightly perturbed by a random 
noise of amplitude e. By this amplitude we specifically mean that a set of uncorrelated Gaussian numbers of zero 
mean and variance was generated, one number for each collocation point in the numerical lattice. Only results for 

> 0 are shown here. The behavior for < 0 is completely symmetrical. 

The initial evolution is well described by the linear stability analysis around the traveling wave Il3yi4|j^ , ^ . 
Typically, as seen from the evolution of the power spectrum, unstable sidebands initially grow. This growth stops 
when an intense competition among modes close to the initial wave and to the broad sidebands is established. 
Configurations during this early nonlinear regime are similar to the ones that would be called riding turbulence and 
described in Section HI. At long times the system approaches one of several possible dynamical states. In general, 
they can be understood in terms of three of them, which are called basic states. In the next section these final states 
are discussed. When the initial winding number is above a critical value Vc, which depends on ci and C 2 , there is a 
transient period between the early competition and the final state during which the winding number changes. 

In Fig. ^ we show in grey levels the phase (p{x,t) for a given run with parameters ci = 2.1 and C 2 = —0.6. The 
space-time defects appear as dislocations in this representation. In Fig. 0b the winding number has been plotted 
as a function of time. The winding number changes from the initial value Vi = 20 to the Hnal value Vf = 14. The 
discrete jumps in v are due to the integer nature of this quantity, and they are smeared out when averages over 
several realizations are performed. The resemblance with the dynamics of the Eckhaus instability of regular waves 
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is striking. In fact, since the changes in v occur on top of a chaotic wave, the analogy is stronger with the Eckhaus 
instability in the presence of stochastic fluctuations [tl7|j3^1 . In the latter case a local wavenumber independent of 
position cannot be defined because of noise, while for phase turbulent waves the disorder is generated by the system 
dynamics. Nevertheless in both cases the configurations can be characterized by a global wavenumber such as k or 
ly. The analo gy i s also instructive since it can be shown that for the one-dimensional relaxational dynamics 

considered in ||37|-|3^ (which is related to Eq. ( 0 ) with Cl = C 2 = 0) there is no long range order in the system, so that 
there is no proper phase transition in the thermodynamic L —*■ oo limit. Despite this, for large but finite sizes and 
long but finite times, sharp transitions are observed and critical exponents and scaling functions can be consistently 
introduced This example should make clear that even in the case that the PT-DT transition would not exist 
in the thermodynamic limit, its characterization in large finite systems is justified. The development of phase slips 
from PT waves of high enough can be viewed as a kind of Eckhaus-like instability for turbulent waves, whereas the 
usual Eckhaus instability appears for regular waves. This similarity was one of the main motivations for the kind 
of analysis that follows. 



FIG. 2. a) Spatiotemporal evolution of the phase coded in grey levels with time running upwards and x in the 

horizontal direction. The lighter grey correspond — —tt and darker to ip{x,t) = tv. The time interval shown in the 

picture goes from t = 500 to 1000 time units of a total run of 10"*. ci = 2.1, C 2 = —0.60, and the initial condition was a TW 
with Vi = 20 that decayed to Vf = 14. The arrow indicates the time at which v begins to change, b) The complete time 
evolution of the winding number for this initial condition. 

For each point in parameter space and initial winding number considered, we have averaged over 50 independent 
random realizations of the white Gaussian perturbation added to the initial wave. Figs. ^ and show the temporal 
evolution of this average v(t) and its variance a for ci = 2.1 and C 2 = —0.83. Four values of the initial winding 
number (vi = 10,15, 20, 25) are shown. Typically, the curve v(t) decays from Vi to a final winding number Vf. 
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FIG. 3. a) Temporal evolution of ^{t) for four different initial winding numbers Vi = 25 (solid), 20(dotted), 15 (dashed) 10 
(dashed-dotted). Ci = 2.1, C 2 = —0.83 (PT regime), b) Winding number standard deviation cr. 

The variance displays the behavior typical of a decay from a unstable state namely a pronounced maximum 
at the time of fastest variation of v{t). The final value of a gives the dispersion in the final values of the winding 
numbers. Although the behavior shown in Fig. is very similar to the observed in for a stochastic relaxational 
case, the scaling laws found there do not apply here. The main qualitative difference is that in a range of Vi the sign 
of the average final v is here opposite to the initial one. In addition for some of the initial winding numbers (i.e. 
Vi = 20 in Fig. ||) v{t) is not monotonously decaying, showing a small recovery after the fast decrease in v. These 
features are also observed for other values of [ci, C 2 ], so that figure || is typical for [ci, C 2 ] in the PT region of Fig. ||. 
For comparison we show v{t) and its variance in Fig. ^for the point Ci = 1.6 and C 2 = —1.0, in the “bichaos” region. 
The main difference is the existence of fast fluctuations in v and cr. They are related to the characteristic dynamics 
of the bichaos regime: The final state depends on the initial conditions and it can correspond to PT, DT or even 
coexistence of both. In the 50 realizations performed all these possibilities were found. When DT appears, there are 
big fluctuations of the winding number around v = 0 that produce the wiggling on the averaged value. More than 50 
realizations should be performed to smooth out such big fluctuations. 
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FIG. 4. a) Temporal evolution of iy{t) for an initial winding number of i/i — 4 in the bichaos regime. ci = 1.6, C 2 = —1.0. 
b) Winding number standard deviation a. 

Returning to the PT parameter regime (Fig. the decay of the initial state is seen to take place during a 
characteristic time which depends on Vi. We quantify this time t as the time for which half of the jump in v is 
attained, r increases as Vi decreases, and there is a critical value of Vc, such that no decay is observed for 
Vi < Vc- Then r diverges (critical slowing down) when Vi approaches Vc from above. This gives a sensible procedure 
to determine v^. Figs. ^ and ||b show l/r as a function of Vi. In Fig. ||a, ci is fixed and the different symbols 
correspond to different values of C 2 . In Fig. C 2 is fixed and the symbols correspond to different values of ci. The 
values of have been estimated by extrapolating to l/r = 0 a linear fit to the points of smallest Vi in each sequence. 
Motivated by we have tried to fit the divergence of r with nontrivial critical exponents, but we have found no 
significant improvement over the simpler linear fit. The values of so obtained are plotted in the insets of Figs. ^ 
and |b. The range of conserved winding numbers [—r'c, r'c] is analogous to the Eckhaus range of stable wavenumbers 
when working below the BFN line. v,. can also be obtained by directly determining the value of Vi below which v{t) 
does not change in any of the realizations. This method can only give integer values of whereas the method based 
on r gives a real number which is preferable when looking for continuous dependences of on system parameters. 
The two methods however give consistent results within the discretization indeterminacy. 




FIG. 5. a) Inverse of the characteristic time for winding number relaxation as a function of the initial winding number. The 
value of Cl is fixed (ci = 2.1) and C 2 varies from near the BFN line (c 2 — —0.48) to the L\ line (c 2 ~ —0.9). Different symbols 
correspond to C 2 = —0.6 (+), C 2 = —0.7 (*), C 2 = —0.75 (O), C 2 = —0.8 (A), C 2 = —0.83 (□). The inset shows the critical 
winding number (vc) as a function of C 2 . b) Idem but the value of C 2 is fixed (c 2 = —0.83) and ci varies from near the BFN 
line (ci ~ 1.33) to ci = 2.5 . Different symbols correspond to ci = 1.6 (+), Ci = 1.8 (*), Ci = 1.96 (O), ci = 2.1 (A), ci = 2.3 
(□), Cl = 2.5 (x). The inset shows the critical winding number {vc) as a function of ci 

The insets of Fig. ||a and Fig. ^ indicate a clear decrease in as the values of Ci and C 2 approach the Li line. In 
fact we know that Vc should be zero to the left of Li, since no wave maintains its winding number constant there. This 
lead us to a sensible method for determining the position of line Li |^, alternative to the one based in the density 
of defects used in 0. It consists in extrapolating the behavior of Vc to = 0. A simple linear fit has been used. The 
same method to determine the line Li has been independently introduced in mi- The coefficients of the linear fit 
are not universal: they depend of the particular path by which the line Li is approached. With this method the line 
Li is determined as the line at which the range of conserved winding numbers [—shrinks to zero. The analogy 
with the Eckhaus instability of regular waves is again remarkable: in the same way as the range of Eckhaus-stable 
wavenumbers shrinks to zero when approaching the BFN line from below, the allowed u range shrinks to zero when 
approaching the Li line from the right. The difference is that below the BFN line the values of the wavenumber 
characterizes plane-wave attractors, whereas above that line, v characterizes phase-turbulent waves. In this picture, 
the transition line PT-DT appears as the BFN line associated to an Eckhaus-like instability for phase turbulent waves. 
For the case of Fig. 5a the PT-DT transition is located at ci = 2.1, C 2 = —0.89±0.02, and ci = 2.60±0.02, C 2 = —0.83 
for the case of Fig. ||b. The agreement with the position of the line as determined by [^j^, where system sizes similar 
to ours are used, is good. For example for ci = 2.1 their value for Li is C 2 = —0.92. The points marked as -I- in Fig. 
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correspond to runs used to determine the position of the transition line Li directly as the line at which defects appear 
in a long run even with i/ = 0. All these ways of determining Li give consistent results. Below the point P, Vc goes 
to zero when the parameters approach the line L 3 , not Li, thus confirming the known behavior that below point P 
in Fig. 1 the line separating phase turbulence from defect turbulence when coming from the PT side is actually L 3 . 

The use of a linear fit to locate the line Li is questionable and more complex fits have been tested. However, the 
simplest linear fit has been found of enough quality for most of the the situations checked. Clearly some theoretical 
guide is needed to suggest alternative functional forms for r'c(ci,C 2 ). We notice that the analogous quantity below 
the BFN line, the Eckhaus wavenumber limit, behaves as ~ for small e, being e the difference between either 
Cl or C 2 and its value at the BFN line. From the insets in Figs, ga or |b, this functional form is clearly less adequate 
than the linear fit used. 

Another interesting point to study is the dependence of the final average winding number on the initial one 
I'i. Fig. I shows an example using ci = 2.1 and C 2 = —0.8. The behavior for other values of the parameters is 
qualitatively similar. Df remains equal to the initial value if r'i < 5 during the whole simulation time, so that Vc ~ 5, 
a value consistent with the one obtained from the divergence of r and plotted in the inset of Fig. ^a. For > Vc, the 
final winding number is always smaller than the initial one. By increasing Vi a minimum on Pj is always observed, 
and then Pf tends to a constant value. Figure]^ also shows the winding number associated with one of the two Fourier 
modes of fastest growth obtained from the linear stability analysis of the initial traveling wave. The one shown is 
the lowest, the other one starts at Pf = 28 and grows further up. Obviously they do not determine the final state in 
a direct way. This is consistent with the observation mentioned above that the winding number instability does not 
develop directly from the linear instability of the traveling wave, but from a later nonlinear competition regime. 



FIG. 6. The final averaged winding number (pf) as a function of the initial one Vi. The initial condition is a TW with 
winding number Ui for ci = 2.1 and C 2 = —0.8. The dashed line corresponds to the lowest of the two Fourier modes of fastest 
growth in the linear regime as a function of Vi. 

As stated in the introduction, a powerful way of distinguishing true phase transitions from effective ones is the anal¬ 
ysis of finite-size scaling We have tried to analyze size-effects from the point of view of v as an order parameter. 
In the DT state such kind of analysis was performed in . Egolf showed that the distribution of the values taken by 
the ever-changing winding number is a Gaussian function of width proportional to '/L. This is exactly the expected 
behavior for order parameters in disordered phases. In the thermodynamic limit the intensive version of the order 
parameter, p/L, would tend to zero so that the disordered DT phase in the thermodynamic limit is characterized by 
a vanishing intensive order parameter. For the PT states to be true distinct phases, the existence of a nonvanishing 
Pc such that p is constant for \p\ < Pc is not enough. The range of stable winding numbers should also grow at least 
linearly with L for this range to have any macroscopic significance. The analysis of the growth of Pc with system size 
has been performed in points ci = 2.1, C 2 = —0.8 and ci = 1.96, C 2 = —0.83 of parameter space, z/c, determined as 
explained before, is plotted in Fig. ^ for several system sizes for which the statistical sample of 50 runs was collected 
for each Pi. 
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FIG. 7. The critical winding number (vc) as a function of the length L of the system is shown. Different symbols correspond 
to Cl = 2.1 and C 2 = —0.8 (A), and ci = 1.96, and C 2 = —0.83 (O). The straight lines are linear hts to the two sets of data. 

There is a clear increasing, close to linear, of as a function of L, thus indicating that for this range of system 
sizes the range of allowed winding numbers is an extensive quantity and then each is a good order parameter for 
classifying well defined PT phases. It should be noted however that for the larger system size for which extensive 
statistics was collected {L = 2048) data seem to show a tendency towards saturation. Thus our study should be 
considered as not conclusive, and larger systems sizes need to be considered. 


III. DIFFERENT ASYMPTOTIC STATES IN THE PT REGION 


Typical conhgurations of the PT state of zero winding number consist of pulses in the modulus |A|, acting as phase 
sinks, that travel and collide rather irregularly on top of the k = 0 unstable background wave (that is, a uniform 
oscillation) SHI- The phase of these conhgurations strongly resembles solutions of the Kuramoto-Shivashinsky (KS) 
equation. Quantitative agreement has been found between the phase of the z/ = 0 PT states of the CGLE and solutions 
of the KS equation near the BFN line fl^. 

For states with zz 0 a typical state n is the one in which an average speed (in a direction determined by the 
sign of v) is added to the irregular motion of the pulses. We have found that in addition to these conhgurations there 
are other attractors in the PT region of parameters. We have identihed three basic types of asymptotic states for 
z/ 7 ^ 0, which we describe below. Other states can be described in terms of these basic ones. Except when exp licitly 
stated, all the conhgurations described in this Section have been obtained by running for long times Eq. (O) with 
the initial conditions described before, that is small random Gaussian noise added to an unstable traveling wave. The 
winding number of these hnal states is constant and is reached after a transient period in which the winding number 
might have chan ged . 

Figs. IJ and pd] show examples of the basic states that we call riding PT (Fig. ||), quasiperiodic states (Fig. ||) 
and frozen turbulence (Fig. [Tl| ). For each hgure: Panel (a) corresponds to a grey scale space-time plot of dx(f(x,t). 
Panel (b) shows the value of this quantity and the modulus of the held (| A |) as a function of position at the time 
indicated by an arrow in panels (a) and (d). Panel (c) shows the spatial power spectrum S(q,t) of A(x,t) for the 
same time. Finally, panel (d) shows the quantity W = J \ dtS{q,t) \ dq, which is a global measure of the temporal 
change in the spatial power spectrum. 

Riding PT. This state (see Fig. ^ is the most familiar one |12 : wiggling pulses in the gradient of the phase 
with a systematic drift in a direction determined by zz. The modulus of the held consists of a disordered spatial 
sequence of small pulses and shocks, with A(x,t) always far from zero. The spatial power spectrum S(q) has a peak 
corresponding to the global wave number k (associated in this case with z/ = —1, so that k = 2'kv/L = —0.012) and 
a broad background associated with the turbulent motion “riding” on the traveling wave. The time evolution of W 
shows a decay towards a huctuating non-zero value, indicating that the power spectrum is continuously changing in 
time as corresponds to the turbulent state reached by the system. 
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FIG. 8. a) Spatiotemporal evolution of dx^pix, t). The lighter grey correspond to the maximum value of dx'^{x, t) and darker 
to the minimum. Last 2000 time units of a run 10"^ time units long for a riding PT state at ci = 2.1 , and C 2 = —0.83. The 
initial condition was a TW with Vi = 20 that decayed to Uf = —1 after a short time, b) A snapshot of | A(x,t) \ and dx^pixA) 
as a function of x for t = 8900 which is indicated by an arrow in a) and d). c) Spatial power spectrum S{q) as a function of 
wavenumber at the same time t = 8900. d)The time evolution of the quantity W defined in the text. The dashed line indicates 
the initial time for picture a). 

Quasiperiodic states. These states (an example is shown in Fig. ^ can be described as the motion of equidistant 
pulses in the gradient of the phase that travel at constant speed on top of the background wave. The fact that the 
periodicity of the pulses and that of the supporting wave are not the same justify the name of quasiperiodic. We show 
later that these states correspond to the ones described in Ref. ||^. In Fig. [^a, the modulus | A | and the gradient 
of the phase clearly exhibit uniformly traveling pulses. The spatial power spectrum S(q) (Fig. W) clearly shows the 
quasiperiodic nature of this state: a central peak, corresponding to the dominant traveling wave, with equally spaced 
peaks surrounding it, showing the periodicity of the pulses. The peaks are not sharp because this configuration has 
been obtained from a random perturbation. The decrease of W in Fig. ^ indicates that the peaks are narrowing. Its 
asymptotic approach to zero indicates that the amplitudes of the main modes reach a steady value and S{q) becomes 
time independent. 
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d) 


FIG. 9. Same as in Fig. ^ but for last 35000 time units of a run 10® time units long for a quasiperiodic state. The initial 
condition is random noise with an amplitude of 0.05. ci = 2.0 , and C 2 = —0.8. b) and c) correspond to a time t = 8 x 10^. 

More perfect quasiperiodic configurations can be obtained from initial configurations that are already quasiperiodic. 
Figure shows the quantity W for a state generated at ci = 2.1 and C 2 = —0.6 from an initial traveling wave with a 
sinusoidal perturbation. The initial traveling wave had Vi = 18 and the winding number of the sinusoidal perturbation 
was v = 22. The travelling wave decayed to a state with Vf = 10 of the quasiperiodic type, cleaner than before. The 
spatial power spectrum (shown in the inset at the time indicated by an arrow in the main picture) shows the typical 
characteristics of a quasiperiodic state. 
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FIG. 10. a)The time evolution of W for a quasiperiodic state. The initial condition is a TW sinusoidally perturbed for 
Cl = 2.1 and C 2 = —0.6 . In the inset the spatial power spectrum S(q) as a function of wavenumber at the time t = 8900 
indicated by an arrow in the main picture. 

Frozen turbulence. This state (see Fig. 0 ) was first described in [^. It consists of pulses in traveling at 
constant speed on a traveling wave background, as in the quasiperiodic case, but now the pulses are not equidistant 
from each other (see Fig. |T^). The power spectrum at a given time is quite different from the one of a quasiperiodic 
state. It is similar, instead, to the power spectrum obtained in the riding PT state: S{q) is a broad spectrum in the 
sense that the inverse of the width, which gives a measure of the correlation length, is small compared with the system 
size. Here however W relaxes to zero, so that the power spectrum finally stops changing (thereby the name frozen). 
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FIG. 11. Same as in Fig. but for last 2500 time units of a run 10'* time units long for a frozen turbulence state. The 
initial condition was a TW of Vi = 12 that decayed to Vf after a short time. ci = 1.75 and C 2 = —0.8. The time of b) and 
c) is t = 8900, indicated by an arrow as in previous hgures. 

This behavior is an indicator of the fact obvious from Fig. |^, that the pattern approaches a state of rigid 
motion for the modulation in modulus and gradient of the phase of the unstable background plane wave. That is, the 
field A{x, t) is of the form: 


A(x, t) = g{x - yt)e^ika:-u.,t+a{t)) 

where g is a uniformly translating complex modulation factor. It is easy to see that configurations of the form S) 
have a time-independent spatial power spectrum. Torcini pl| ] noticed in addition that the function a{t) is linear in t 
so that the solutions are in fact of the form 

A(x, t) = f{x - (3.2) 

where again f(x — vt) is a complex valued function and lo can differ from / and g differ only in a constant phase. 
The envelopes g{x — vt) or f{x — vt) turn out to be rather irregular functions in the present frozen turbulence case, 
whereas they are periodic in the quasiperiodic conhgurations discussed above. 

After presenting the basic states, we continue addressing some interesting mixed states that can be described in 
terms of them. Most of the configurations ending up in the frozen turbulence or in the quasiperiodic states have long 
time transients of the riding turbulence type. Only at long times a decay to rigid propagation occurs. There are cases 
in which a different type of decay happens. For example Fig. ^ shows a case in which the system jumps from a very 
strong riding turbulence regime to another state, also of the riding turbulence type, but much more regular. The 
quantity W, shown in Fig. |T^, turns out to be a valuable tool in distinguishing the different regimes: a superficial 
look at Fig. |l^ could be easily misunderstood as indicating the approach of the system towards a frozen turbulence 
state, but the lack of decay towards zero of W identifies the final state as riding turbulence. The arrows indicate the 
jump to the second state. Fig. |l^ shows a state characterized by a recurrence between two different riding turbulence 
states, showing a kind of temporal intermittency. 
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FIG. 12. a) Spatiotemporal evolution of dx^{x, t) for a riding turbulence state that decays onto another one. Ci = 2.5, 
C 2 = —0.75. The initial condition was a TW of Vi = 20 that decayed to Vf = —2 in a short time, b) Time evolution of W. 
The dashed lines indicate the time interval shown in a) (from t\ = 2500 to t 2 = 6500 of a run 10® time units long). The arrow 
indicates the transition from one of the riding turbulence regimes to the other one. 



0 aooo 4000 6000 BOOOIOOOO 


FIG. 13. a) Spatiotemporal evolution of showing intermittency between riding turbulence states. ci = 2.1, 

C 2 = —0.83. The initial condition is a TW of Vi = 1 that did not change. b)Time evolution of W. The dashed lines indicate 
the time interval shown in a) (from ti = 1000 to t 2 = 8500 of a run 10'^ time units long). The arrows indicate the end of a 
riding turbulence regime and the beginning of another one. 


Finally Fig. shows a riding turbulence state with zero winding number. This is not however a typical configu¬ 
ration, since usually for zz = 0 there is no preferred direction for the pulses to drift, whereas the figure shows that in 
fact there is a local drift at some places of the system. It turns out that this state can be understood as composed by 
two domains of different local winding number: v = 1 and v = —1, so that globally v = 0. The pulses travel either in 
one direction or in the other depending of the region of the system in which they are. In Fig. & a snapshot of the 
gradient of the phase dx^{x, t) and the phase itself ^p{x^ t) is shown. Lines showing the average trend in the phase are 
plotted over the phase, clearly identifying the two regions in the system. This coexistence of the different basic states 
at different places of space, or at different times as in Fig. 13, was already mentioned in ||ll|] where it was argued to 
give rise to a kind of spatio-temporal intermittent behavior. 
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FIG. 14. a) Spatiotemporal evolution of dx^{x,t). The time interval corresponds to 5500 to 7500 time units of a run 10^ 
time units long for a riding PT state at ci =2.1 and C 2 = —0.83. The initial condition was a TW with Ui = Q that did not 
change, b) Snapshots of ip{x,t) and dxg>{x,t) as a function of x at the time t = 6980 indicated by an arrow in a) and d). The 
dashed lines in the graph of g>{x,t) indicate average slopes, that is “local” wavenumbers, c) Spatial power spectrum S{q) as a 
function of wavenumber at the same time t = 6980. d)Time evolution of W. Dashed lines indicate the time interval of picture 
a). 


Given the large variety of configurations that are observed, and the very long transients before a jump from one 
state to another occurs, it would be difficult to conclude from numerical evidence alone that the three kinds of states 
considered as basic above are true asymptotic states. Some analytical insight would be desirable to be sure that these 
three states are attractors of the dynamics. The next Section is devoted to provide such analytical justification. 


IV. ASYMPTOTIC STATES IN TERMS OF THE PHASE DYNAMICS 


The question on whether it is possible or not to describe the PT regime of the CGLE from a closed equation for 
the phase alone has been posed by several authors 16,3^, ^ , ^ . A phase equation is obtained by considering a 
long wavelength perturbation of a plane-wave solution in the GGLE (|l . l| ). It is clear that this phase equation will 
only describe phase dynamics close to the homogeneous plane-wave (that is the one with z/ = 0) if the perturbation is 
made around the spatially-homogeneous solution. In order to get a description of PT at 7^ 0 the expansion should 
be done for a perturbation on a traveling wave solution with wavenumber (fc) different from zero, 


A = (x/l-P + a(x, , (4.1) 

Here k is taken as k = If A satisfies periodic boundary conditions the same conditions apply to (f> because any 
global phase winding is included in k (the total phase is ip = kx + cj)). From general symmetry arguments the general 
phase equation for fc yf 0 should read, up to fourth order in gradients: 

dt4> = flo - Vgdx4> - D2dl(l) + + D^dlcj) -b Di2{dx(j)){dl(j)) 

- + D^^{dx^){dlrj>) A D22{dlrj>f + D^i2{dx^f{dlcfy P ... . (4.2) 
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When Vg = D 3 = D 12 = D 13 = D 22 = D 112 = 0, Eq. ( |4.2| ) reduces to the Kuramoto-Sivashinsky (KS) equation 
that is the lowest order nonlinear phase equation for the case fc = 0. For k ^ 0, Eq. (O) was syste mati cally 
derived up to third order in gradients in |^. An easy way of obtaining the values of all the coefficients in (4.2) was 
discussed in M: First, Qq is related to the frequency of the plane-wave solutions: 


Oo = -Wfc = -C 2 - (ci - C 2 )k^ 


(4.3) 


Second, the linear terms can be obtained from the eigenvalue X(k, q) corresponding to the phase-like branch in the 
linear stability analysis of the wave of wavenumber k with respect to perturbations of wavenumber q : 

A(fc, q) = -ivgq + D 2 q^ - W^q^ - D^q'^ + 0{q^) (4.4) 


with 


D 2 

D 3 


Di 


2k{ci — C 2 ) 

\ 2k'^{l + cl) 

2fc(l -I- C 2 ) [—Cl -f (ci -I- 2 c2)A:^] 


1 


(i-k^y 

{cf (1 -f cl) - 2k'^{l + cl){cl + 6 C 1 C 2 ) 


2(1- fc2)3 

k^ [4 -f (1 -f cl){cl + 12 ciC 2 ) -f c^(24 -f 20c^)] } . 


(4.5) 

(4.6) 

(4.7) 


(4.8) 


Third, the nonlinear terms can be obtained from the following consistency relationship: If ( [4.lD is an exact solution 
of the CGLE, then (j){xy) satisfies the phase equation with coefficients depending on k. In addition if (-^1 — 4 - 

ai(a;, jg another exact solution of the CGLE, then satisfies a similar phase equation but with 

coefficients depending on ki instead of k. But this solution can be written as (-^/l — -|-ai(a:, i))edfe 2 :-i-(fci-fe)x+ 0 i(x,t)) 

so that {ki — k)x + (j)i{x,t) is also solution of the phase equation with coefficients depending on k (with different 
boundary conditions). By combining the two equations satisfied by (f>i and expanding the coefficients depending on 
ki as a power series around k (assuming ki — k small) the following relationships between linear and nonlinear terms 
are obtained: 


So that 


-Dll = - 


^ dVg 

2 dk 


Di 2 


dD 2 n n 


(4.9) 


-Dll — C2 — Cl 


-D 12 = 
-D 13 = 
D 112 = 


4fc(l -|- cl) 

~ (1 -P)" 

y^l-o^+<io,k■‘ + (2c,+c,)k‘] 

[l — k^) 

2(l + ci)(3P + l) 

( 1 - fc 2 )^ 


(4.10) 

(4.11) 

(4.12) 

(4.13) 


The coefficient D 22 is only obtained following the method to higher order in (fci — fc). The coefficients up to third 
order in gradients can be found also in and approximate expressions for them are given in . 

The traveling wave of wavenumber fc be com es unstable when the coefficient D 2 becomes positive. One expects that 
the first terms in the gradient expansion (H) give a good description of the phase dynamics in the weakly nonlinear 
regime, that is D 2 positive but small (note that for a given fc 7 ^ 0 this includes part of the region below the BEN line 
in Fig 1). The arguments presented in imply that the relative importance of the different terms in a multiple scale 
expansion in which D 2 is the small parameter can be established by considering (j) ^ dx ^ ■ Then the dominant 

terms close to the instability of wave fc are the ones conta inin g Dq and Vg. After them, the terms with coefficients 
D 3 and £>11 are the most relevant. Up to this order Eq. (^ ) is a Korteweg-de Vries equation (KdV). The terms 
with D 2 , -D 4 and D 12 appear at the next order. The importance of the terms in D 2 and D 4 for a qualitatively correct 
description of phase dynamics is obvious since they control the stability properties of the wave of wavenumber fc. The 
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importance of the term with coefficient D 12 was stressed in 00 : if it is large enough it can change the character 
of the bifurcation from supercritical to subcritical. 


The detailed comparison of the reduced dynamics (4.2) with the complete CGLE phase dynamics is beyond the scope 
of the present paper. The aim of this Section is to use Eq.(4^) just to get some understanding of the asymptotic 
states presented in Section III. To this end we will use the detailed results available from the work of Cha ng et 
al. Q. These results are obtained for the so-called Kawahara equation which is Eq. (|4.2| ) with 

Di 2 = Di 3 — D 22 = D 112 = 0. The term D 12 , which according to Kuramoto estimations |^| is of the same order 
for small D 2 as the terms in D 2 and £> 4 , will thus be neglected. It would be certainly necessary to consider the 
modifications introduced by the term D 12 into the results of [Q. This will be briefly discussed at the end of this 
section. At this point it is interesting to note that, to our knowledge, the only quantitative comparison of the phase 
dynamics with k ^ 0 obtained from a phase equation and from CGLE is E&- But the phase equation used in 
these references is the one presented in |^, in which the nonlinear terms considered are only those with coefficients 
Dll and D 13 . In addition Dii,Di 3 , and the coefficients of the linear terms are considered only up to first order in k. 
Despite these limitations, in particular the absence of the D 12 term, the phase equation is found to reproduce well the 
phase dynamics of the CGLE, an agreement that degrades when the term in £ >13 is suppressed [0|. Clearly further 
work is needed to establish firmly the relevance of the different terms in (4.2) [^. Our study will be restricted to 
the situation of 0 | (that is £>12 = £>22 = D 13 = D 112 = 0 ) since no study of comparable detail for a more complete 
equation is available in the literature. 

The situation of interest here is the one in which the traveling waves are unstable against a finite band of wavenum¬ 
bers, so that £> 2 , D 4 > 0. Making the following changes of variables in (4.2) with D 12 = D 22 = D 13 = £*112 = 0: 



DuDI^^ 
HX.M =- 

2 dI'^ 


the Kawahara equation 50 Eiliil is obtained 


drU = — 4ud^u — 5d^u — d*u 


with 


5 = - 


D. 


D2 D4 


(4.14) 


(4.15) 


(4.16) 


Since </> is periodic in x, u{x,t) is periodic in y. In addition u(x,T)dx = 0. To have some intuition on the meaning 
of the parameter S, its expansion at small k reads 


S « 2 V 2 k sign(ci)i 


1 +cl 


1 -I- C 1 C 2 


C>(P) . 


(4.17) 


It should b e not ed that 6 does not diverge at the BEN line, as the expansion (4.17) seems to suggest, but below it. 
Erom Eq. (4.16) it is clear that S diverges where D 2 vanishes indicating that the corresponding traveling wave of 
wavenumber k has become Eckhaus unstable. 


The Kawahara equatio n (4.15|) has been considered in the context of surface waves on fluid Aims falling down an 
inclined or vertical plane |54 ], and also as a simple renerafization of the KS or the KdV equations [ 0j5^ . It has also 
been considered in the context of growth shapes 1^. It reduces to KS for <5 = 0 (or equivalently for k = 0) when 
written in th e ori ginal variable (p. 

Equation ( 4.15| ) has periodic, soliton-like, spatially-irregular, and spatio-temporally chaotic solutions. 0-0. In 
fact, all of these solutions have been analytically shown to exist 0]. All of them except the isolated soliton-like 
solution 01 are stable in some parameter regimes |0j . These kinds of solutions should manifest themselves (provided 
the approximate phase description holds) in the time evolution of the phase gradient {= k + d(j)) of the solutions 
of the CGLE 0 in the PT regime. The analytical results in thus provide a Arm basis for true existence of the 
numerically observed states described in Section IB. 

The detailed bifurcation analysis in 0] also gives detailed predictions for the wound states of the CGLE, within 
the range of validity of the phase description. We will reproduce here some of the results in [Esl and reinterpret them 


16 


























in terms of the gradient of the phase of CGLE solutions. Our interest is centered in the rigidly moving train of pulses 
( frozen turbulence an d qu asiperiodic states ) obs erved in several of the numerical simulations reported in Section 01. 
They are of the form (|3.2D, and because of (4.14) we have 


i(x,'r) = H{^) , 


(4.18) 


with ^ = X ~ be ing v the velocity of the train of pulses we want to describe in units of x r. The partial 
differential equation ( 4.15| ) is reduced to an ordinary differential equation (ODE) for 

(4.19) 


The primes denote differentiation with respect to After an integration: 

H'" + 5H" + H' -vH + 2H^ = Q . 


(4.20) 


Q is fixed in a nontrivial way by the condition J Hd^ = 0 which follows from our periodic boundary conditions. This 
third order ODE can be rewritten as a three-dimensional dynamical system: 


Ui = U2 
/ 


U2 = W3 

U 3 = CUi — U2 — Sus — 2(ui)^ 


(4.21) 


with 


M0 = h{x)-1 + \I^ + ^, 


c = vSQ + zA . 


(4.22) 


Different qualitativ e beh aviors in phase space of the solutions of the dynamical system (4.21) are related to the shape 
of the solutions of ( 4.1£ ) |Q. This is illustrated in Fig. 
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FIG. 15. Schematic relationship between traj ector ies of the dynamical system (1.21) in phas e space (left column) and 
solutions u(x,t) = = X ~ of equation (4.15) (right column), al) Fixed point of ( [4.2l| ), a2) uniform solution of 

(4.15) (traveling wave in the CGLE (|l.l|) ). bl) Periodic solution (limit cycl e) of (4.21), b2) per iodic train solution of (4.15) 
(qu asipe riodic solution of the CGLE ([l.ip). c l) H omoclinic trajectory of (4.21), c2) single pulse of (4.15). dl) Chaotic trajectory 
of (4.21), d2) spatially irregular solution of (4.15) (frozen turbulence in the CGLE). 


We stress that all the solutions o f ([4.19| ) represent uniformly translating solutions of ( 4.15| ). No information is given 
on more complicated solutions of (4.15). The left column of Fig. O shows the possible trajectories of the dynamical 


system (4.21) while t he rig ht column shows the co rresp onding solution of (4.2C), or equivalently tt(Y, t) = g (£ = 
Y — vt ) in equ ation ( 4.15 ). For a fixed point in ( 4.2l|) (fig. |l5|al ) we get a homogeneous soluti on in ( |4.15|) ( fig. 
p^2) and (via (4.14)) a traveling wave solution in t he C GLE ( |l.lD . For a periodic trajectory in ( |4.2lD (fig. |l5|bl ) 
we get a train of periodic p ulses in the solution of ( 4.15 ) (fig. 1l^2) and a quasiperi odic s olution in CGLE (^J). 
An homoclinic trajector y in ( p3I| ) (fig . |l5|cl) corresponds to a single pulse solution in ( 4.15| ) (fig. p^c2). Finally for 
a chaotic trajectory in ( 4.21 ) (fig, [l^ l) we have an irregular solution H {^) th at corresponds to a rigidly traveling 
spatially irregular solution of ( 4.15| ) (fig. p^2). The chaotic solutions of (4.21) are of the Shil’nikov type |Q. This 
means that the disordered configurations (and thus u and dx^) consist on nearly identical pulses irregularly 

spaced. This corresponds to the state named frozen turbulence for solutions of the CGLE. 

The detailed analysis of is done on the one hand by following the sequence of bifurcations of t he sta te in which 
iL is a constant and of the state in which H is close to the KdV soliton (with adequate rescaling Eq. (4.15) reduces to 
the KdV in the limit 6 —> oo). On the other hand the powerful glob al theorems of Shil’nikov and their generalizations 
[^8|-|^ are used to establish the structure of the solutions of ( 4.21 ). The results of Q relevant to our purposes can 
be summarized as follows (they can be read-off from figure 3 of Ref. ): 


1. Periodic solutions of ( |4.2l| ) exist for all values of 6 provided | c |>| 5 |. They are organized in a variety of 
branches. Solutions in the same branch differ by their periodicity, and each branch ends in a different kind of 
solitary-wave solution (infinite spatial period). The shape of the different solitary wave solutions characterizes 
the different branches. 

2. For I 5 I ^ 1.1 only one of the branches of periodic solutions (the main branch) remains. 


3. Chaotic solutions to (4.21) exist only for | S 0.84. 


In addition Chang et al. |Q obtained results also for the full equation ( 4.15| ), without the restriction to rigidly 
traveling waves. Their numerical and analytical results can be summarized as 


4. Periodic solutions in the main branch with its wavenumber within a given range are linearly stable for all S. A 
more precise determination of the range of stable wavenumbers for large 6 was performed in . 

5. In addition to the periodic solutions there are also spatio-temporal chaotic attractors for all S. 

6 . If I 5 |> 1.1 only two of these strange attractors remain. For | J |> 3 their basin of attraction seems to be much 
smaller than the one of the periodic solutions. 


Expression ( 4.16 ) with (4.5)-( 4.13 ) gives the relation between S and the parameters of the CGLE. | 5 |= oo 
corresponds in Fig. 1 to the line at which the wave of wavenumber k becomes Eckhaus unstable. It is approximately 
parallel and below the BFN line. The other lines of constant 6 , for fixed k, are also approximately parallel to the BFN 
line, and decreasing | 6 \ corresponds to entering into the PT region and going deep into it. All these lines concentrate 
onto the BFN line as k approaches zero: for k = 0, S = 0 except on the BFN line 1 -I- ciC 2 = 0 where <5 is undefined. 
We now rephrase the conclusions above in terms of the th ree basic asymptotic states of the CGLE in the PT regime. 
They will be valid as long as the phase description (4.15) remains accurate. 


1. There are PT solutions of the quasiperiodic type for all values of the parameters (as long as the phase description 
remains valid). Bounds on their velocity can be in principle obtained, but this is nontrivial since Q is only known 
in an implicit way. 


2. Increasing | <5 | by approaching the Eckhaus instability for a given k {D 2 = 0), or by increasing the winding 
number reduces the variety of quasiperiodic solutions. 

3. Frozen turbulence solutions exist only for | S 0.84, that is far enough from the line D 2 = 0 or for small enough 
winding number. 


4. There are linearly stable solutions in the main quasiperiodic branch for all values of the parameters. 
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5. There are also riding turbulence attractors for all values of parameters. 

6 . For I 5 |> 3, that is at high winding number or close enough to the line 132 = 0 the quasiperiodic solutions have 
a basin of attraction larger than the riding turbulence ones. 

A general feature of these conclusions is that the important quantity is D 2 , that is the distance in parameter space 
from the line at which the fc-wave became Eckhaus unstable. This line is below the BEN line for k ^ 0. Thus not 
only traveling waves, but also quasiperiodic, frozen turbulence, and riding turbulence attractors should exist below 
the BEN line for k ^ 0. In practice it is relatively easy to find quasiperiodic states below but close the BEN line, 
but we have been unable to find the other two states so far. The difficulty in finding riding turbulence states can be 
a consequence of the small range of winding numbers for which they are stable {\i^\ — L\k\/{2Tr) < Vc) so that the 
observability condition | (5 |< 3 immediately brings us above the BEN line. Another possibility is that the instability 
of the V = 0 plane-wave attractor at the BEN line has consequences of a global character beyond the validity of the 
phase description. 

The above predictions imply that the more promising zone for obtaining quasiperiodic solutions starting from 
random perturbations on a traveling wave of given winding number is for parameter values close and above D 2 = 0, 
or for high winding number (| i5 |> 3). In any case no frozen turbulence should be observed in that zone. 

Some qualitative aspects of the conclusions above have been shown to be correct. In particular Torcini and col¬ 
laborators 0,11 have shown that the average maximal Lyapunov exponent, quantifying the proportion of initial 
conditions that fall into the spatio-temporal chaotic strange attractors, is a decreasing function of v. 

Our numerical solutions also agree with the prediction that quasiperiodic solutions show up more easily for small 
D 2 . However, their basin of attraction appears to be much smaller than the implied by the conclusions of the phase 
description since it is reached with very low probability from our initial conditions. This is specially true above the 
BEN line. The reason for this is probably the effect of the neglected term D 12 , which is known to reduce the range 
of stable periodic solutions and even to eliminate it by making the bifurcation subcritical ||^. Above the BEN 
line the attractor that we observe more frequently at high winding number from our initial conditions is the frozen 
turbulent state. 

A more detailed checking of the predictions above would be desirable. This is however beyond the scope of the 
present paper since a detailed theoretical analysis of the global properties of the phase space for the equation containing 
the term D 12 wo uld be probably needed beforehand. A promising alternative can be the study of the exact equation 
for f{x — vt) in ( |3.2| ) obtained in [^. 


V. FINAL REMARKS 


One-dimensional wound-up phase turbulence has been shown to be much richer than the case v = Q. The main 
results reported here, that is the existence of winding number instability for phase-turbulent waves, the identification of 
the transition PT-DT with the vanishing of the range of stable winding numbers, and the coexistence of different kinds 
of PT attractors should in principle be observed in systems for which PT and DT regimes above a Hopf bifurcation 
are known to exist |Q. To our knowledge, there are so far no observations of the ordered PT states described above. 
There are however experimental observations of what seems to be an Eckhaus-like instability for non-regular waves 
in the printer instability system This suggests that the concept of a turbulent Eckhaus instability can be of 

interest beyond the range of situations described by the CGLE. A point about which our study is inconclusive is the 
question on the existence of PT in the thermodynamic limit. The identification of v as an order parameter identifies 
the continuation of Fig. towards larger system sizes as a way of resolving the question. It should be noted however 
that although a linear scaling of the order parameter with system size is usual in common phase transitions, broken 
ergodicity phase transitions, as the present one, generate usually a number of ordered phases growing exponentially 
with L, not just linearly p3J^. We notice that the results of Section III show that the states of a given v are not 


pure phases, but different attractors are possible for given v. An order parameter more refined than v should be able 
to distinguish between the different attractors and, since some of them are disordered, the result of an exponentially 
large number of phases at large L would be probably recovered. The results presented in Section 0 give a justification 
for the existence of the several wound states observed, an especific predictions have been formulated on the basis of 
previous anal itical and numerical results. Further work is needed however to clarify the importance of the different 
terms in (O) and the validity of a phase description. 
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APPENDIX A: NUMERICAL INTEGRATION SCHEME 


The time evolution of the complex field A{x, t) subjected to periodic boundary conditions is obtained numerically 
from the integration of the CGLE in Eourier space. The method is pseudospectral and second-order accurate in time. 
Each Fourier mode Aq evolves according to: 


dtAq(t) = -aqAq(t) + $,(t) 


(Al) 


where aq is (1 -I- ici)q^ — 1, and ^q is the amplitude of mode q of the non-linear term in the CGLE. At any time, the 
amplitudes are calculated by taking the inverse Eourier transform A{x,t) of Aq, computing the non-linear term in 
real space and then calculating the direct Fourier transform of this term. A standard FFT subroutine is used for this 
purpose M- „ 

Eq. (Al) is integrated numerically in time by using a method similar to the so called two-step method Eor 


convenience in the notation, the time step is defined here so as the time is increased by 2dt at each iteration. 

When a large number of modes q is used, the linear time scales aq can take a wide range of values. A way of 
circumventing this stiffness problem is to treat exactly the linear terms by using the formal solution: 


Aq{t) = e ^Aq{to)e‘^Ao ^ $,(s)e“‘'"ds^ . 


Erom here the following relationship is found: 

Aq(t + St) Aq{t — St) 


nt-\-5t 


^ — OLnSt 




5t 


= e ^ 


$,(s)e“‘>'*ds. 


' t—St 


The Taylor expansion of of 4)q(s) around s = t for small St gives an expression for the r.h.s. of Eq. (A3): 




^OtqSt _ g CXq6t 


O(St^) . 


(A2) 


(A3) 


(A4) 


Substituting this result in {\K^ ) we get: 




1 — e 


-20CaSt 


-<i>q{n) + 0{St^) . 


(AS) 


where expressions of the form /(n) are shortcuts for f(t = nSt). Expression (AS) is the so called ’’slaved leap frog” 
of Frisch et al. |^. To use this scheme the values of the field at the first two time steps are required. Nevertheless, 
this scheme alone is unstable for the CGLE. This is not explicitly stated in the literature and probably a corrective 
algorithm is also applied. Obtaining such correction is straightforward: Eollowing steps similar to the ones before one 
derives the auxiliary expression 


A,(n)=e-“-^*A,(n-l) + 


1 — e 




-^q{n-l) + 0[St^), 


(A6) 


The numerical method we use, which we will r efer to as the two-step method, provides the time evolution of the 
field from a given initial condition by using Eqs. (|^ ) and (^) as follows: 

1. <i)g(n — 1) is calculated from Aq{n — 1) by going to real space. 

2. Eq. ([A^ ) is used to obtain an approximation to Aq{n). 

3. The non-linear term $q(n) is now calculated from this Aq(n) by going to real space. 
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4. The field at step n + 1 is calculated from (A5) by using Aq{n 


1) and 4)g(n). 


At each iteration, we get Aq{n+1) from Aq(n—1), and the time advances by 2dt. Note that the total error is 
despite the error in the intermediate value obtained with Eq. ( [^ ) is 0{St^). The method can be easily made exact 
for plane waves (|l.2| ) of wavenumbe r k (and then more precise for solutions close to this traveling wave) simply by 
replacing the nonlinear term <i>g in & by + (1 + JC 2)(1 — k‘^)Aq, and subtracting the corresponding term from 
aq. We have not implemented this improvement because we were mostly interested in solutions changing its winding 
number, so that they are not close to the same traveling wave all the time. 

The number of Fourier modes depends on the space discretization. We have used dx = 1 and usually N = 512. 
The time step was usually dt = 26t = 0.01. The accuracy of the numerical method has been estimated by integrating 
plane-wave solutions. The amplitude and frequency of the field obtained numerically will differ slightly from the exact 
amplitude and frequency, not only due to round-off errors, but also due to the fact that the method is approximate. 
The method has been tested by using a stationary unstable traveling wave of wave number k as initial condition. 
The numerical errors will eventually move the solution away from the plane-wave unstable state. To be precise, in a 
typical run with ci = —1.0 and C 2 = 2.4, with dt = 0.01 and k = 0.123, the amplitude was kept constant to the fifth 
decimal digit during ~ 8000 iterations. In comparison, when a Gaussian noise with an amplitude as small as 10“^ 
is added to the traveling wave, the modulus is maintained equal to its steady value (up to the fifth decimal) during 
1500 iterations. The frequency ujq determined numerically by using dt = 0.01 fits the exact value up to the fourth 
decimal digit. 

The integration method introduced here has also been applied succesfully to the case of two coupled equations, or 
equivalently a Vectorial CGLE j^,^. 
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